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Abstract. We investigate the implications of a massive 
binary system in the centre of the gamma-ray blazar 
Mkn 501 and show that the periodical behaviour recently 
observed in the TeV and X-ray lightcurves may possibly be 
related to the orbital motion of the relativistic jet emerg- 
ing from the less massive black hole. For the special rel- 
ativistic jet properties inferred from emission models, we 
derive an intrinsic orbital period of (6 — 14) yrs and a 
centre-of-mass distance of (2.0 — 3.5) x 10 16 cm. If the bi- 
nary is very close with a separation of the order of that for 
which gravitational radiation becomes dominant, we find 
a maximum primary mass of ~ 10 8 Mq and a correspond- 
ing secondary mass in the range of ~ (4 — 42) x 10 6 Mq 
depending on the intrinsic jet properties. Such values are 
in line with the black hole masses expected from merger 
scenarios. 

Key words: Galaxies: active - BL Lacertae objects: in- 
dividual: Mkn 501 - Galaxies: jets 



1. Introduction 

Binary black hole systems (BBHSs) are expected to be 
common in the universe as a result of mergers between 
galaxies. In the underlying picture for the morphological 
evolution, galaxies were formed as part of a hierarchical 
clustering process (e.g. White 1997). Giant elliptical galax- 
ies, such as the host galaxy of Mkn 501, appear to be the 
products of mergers between spiral galaxies (cf. Fritze v.- 
Alvensleben 1996). Since the brightest galaxies generally 
seem to contain massive black holes in their nuclei (e.g. 
Rees 1984; Kormendy & Richstone 1995; Ho 1998; Magor- 
rian et al. 1998; Richstone et al. 1998), merging would nat- 
urally lead to the formation of massive BBHS (Begelman 
et al. 1980, abbreviated: BBR 80; Rees 1994; Artymow- 
icz 1998; Richstone 1998). The so formed binary black 
hole is expected to spent most of its time at a separation 
of - 0.1 - 1 pc for masses of 1O 8 M (BBR 80). How- 
ever if the binary looses further angular momentum e.g. 
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by slingshot interaction with new stars from subsequent 
merging events (Roos 1988; Roos et al. 1993), infall of gas 
(BBR 80; cf. also Gould & Rix 2000) or by interactions 
with an accretion disk (Ivanov et al. 1999), gravitational 
radiation will eventually become important and the binary 
evolution could proceed rapidly to coalescence. 

Up to now, several phenomena have been attributed to 
BBHSs: e.g. misalignment (cf. Conway & Wrobel 1995), 
precession (BBR 80) or wiggling of jets, where the latter is 
supposed to be induced by the orbital motion (Kaastra & 
Roos 1992; Roos et al. 1993). Periodic outburst activity in 
the quasar O J 287, has commonly been related to a BBHS 
and is thought to arise due to tidal perturbation (Sillanpaa 
et al. 1988) or due to one black hole crossing the accretion 
disk of the other (Letho & Valtonen 1996). Other BBHS 
scenarios assume a pair of bent jets (Villata et al. 1998) or 
the precession of the disk under the gravitational torque 
(Katz 1997). 

In the particular case of Mkn 501, the complex mor- 
phology of its radio jet and the peculiar behaviour of its 
spectral energy distribution (SED) have prompted elabo- 
rate models relating these properties to a BBHS: Conway 
& Wrobel (1995), for example, have proposed a saturated 
helix model in order to explain the misaligncment of the 
radio jet on parsec and kiloparsec scales. Villata & Rai- 
teri (1999) have argued that the X-ray variations in the 
SED of Mkn 501 might be due solely to the changing ori- 
entation of a helical synchrotron emitting jet in a close 
BBHS. The recent discovery of periodicity in the TeV and 
X-ray fluxes believed to be associated with moving fea- 
tures in the jet of Mkn 501 might add another aspect for 
assessing the relevance of a BBHS in this galaxy. 

Mkn 501 is one of at least four active galactic nuclei 
which have been detected at TeV energies (for review, 
see Catanese & Weekes 1999). Being the second closed 
among these with a redshift of z = 0.034, Mkn 501 has 
been historically classified as an X-ray selected BL Lac 
object showing virtually no emission lines, and is hosted 
by the elliptical galaxy UGC 10599 (Stickel at al. 1993). 
As a BL Lac object, Mkn 501 belongs to the blazar class of 
AGN which are thought to have relativistic jets oriented 
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at a small viewing angle, thus yielding a strong Doppler 
enhancement of the observed flux. 

At the beginning of 1997, Mkn 501 had suddenly un- 
dergone a phase of high activity becoming the brightest 
source in the sky at TeV energies. Subsequent multiwave- 
length campaigns revealed a variable, two component SED 
with a low energy part extending up to 100 keV (Pian et 
al. 1998) and a high energy part which extends at least 
up to 20 TeV (Samuelson et al. 1998; Konopelko 1999). 
During this activity phase, particular types of variability 
have been observed (e.g. Protheroe et al. 1998), consisting 
of flaring episodes of several days and additional intraday- 
variabilities. While the TeV and X-rays variations seem to 
be well correlated, the evidence for correlations with the 
optical U-band appears to be rather weak (e.g. Catanese et 
al 1997; Djannati-Atai et al. 1999; Aharonian et al. 1999). 
One of the most fascinating features is the observed pe- 
riodicity in the TeV region with a period in the range of 
(23 — 26) days, which has been found in the data taken 
by several Cherenkov telescopes (see Protheroe et al. 1998; 
Hayashida et al. 1998). Additionally, a recently performed 
analysis of RXTE-ASM X-ray data from April to July 
1997 also seems to support a periodicity of ~ 23 days 
(Kranich et al. 1999; Nishikawa et al. 1999). 

Here, we speculate on the possibility that this period- 
icity in the flaring state arises due to the orbital motion 
of the relativistic jet in a BBHS where the nonthermal 
radiation is emitted by a relativistic jet which emerges 
from the less massive black hole, the periodicity thus be- 
ing due mainly to geometrical origin (i.e. Doppler-shifted 
modulation) . 

2. On the influence of a BBHS in Mkn 501 

2.1. Doppler-shifted flux modulation 

Consider a simple binary model for Mkn 501, where the 
binary orbit is is assumed to be circular because dynamical 
friction between two parent galaxies during the merger 
might ensure that the initial eccentricity of the resulting 
binary is small (e.g. Polnarev & Rees 1994). Using Kepler's 
third law, the intrinsic orbital frequency fi/. of a binary 
with separation d is given by 



y/G (m + M) 

W 2 



(1) 



where m and M denotes the mass of the smaller and the 
larger hole, respectively, and G is the gravitational con- 
stant. 

Let us further assume that the observed jet is formed 
by the less massive black hole and that the nonthermal 
X-ray and 7-ray radiation in the flaring state is emitted 
by a relativistic emission region (e.g. knot, blob, shock) 
which propagates outwards from the core along the jet 
with gamma factor 7^. Owing to the (non-relativistic) or- 
bital motion, the true trajectory of the knot is a long 



stretched helix. The modulation of the emission then oc- 
curs as a consequence of the slight variation of the in- 
clination angle i due to the orbital 0-component of the 
velocity field of the knot. This observed flux modulation 
by Doppler boosting is well-known and can be written in 
the form 



S(v) =S 3 S'{v')=5 3+a S'{v) 



(2) 



where S' is the spectral flux density measured in the co- 
moving frame, S(t) — l/(7fc[l — Pt cos9(t)]) the Doppler 
factor, 9(t) is the actual angle between the velocity (3b = 
xi>(t)/c of the emission region and the direction of the ob- 
server, and where the final equality holds if the source has 
a spectral index a. 

Due to the orbital motion around the center-of-mass, 
the Doppler factor for the emission region is a periodi- 
cal function of time. In the simplest case where the angle 
between the jet axis and the direction of the total angu- 
lar momentum of the binary is assumed to be zero (e.g. 
neglecting any kind of precessional motion) the Doppler 
factor may be written as 



1 — (v z cos i — QkR sin i sin Qj. t) jc 



(3) 



with R = M d/ (m+M) being the centre-of-mass distance, 
v z the outflow velocity in the direction of the total angu- 
lar momentum, i the inclination between the jet axis and 
the line of sight and c the velocity of light. Obviously, 
the Doppler factor becomes maximal for t = 0.75 Pk and 
minimal for t = 0.25 Pk, where Pk = 2 it /ilk denotes the 
keplerian period. 

From the TeV flux ratio of / ~ 8 between the maxi- 
mum and the minimum state during the observation (cf. 
Protheroe et al. 1998; Hayashida et al. 1998, Aharonian et 
al. 1999) and the assumption that the periodicity arise in 
the main due to geometrical origin, we now obtain the con- 
dition <5 max /<5 min ~ /V(3+«) ( se e Eq. (§)). Consequently, 
by using Eq. (0) one finds 



ilk R 



f 



l/(3+a) 



1 



/l/(3+a) + 1 



1 Vz 

— cot I 

sini c 



(4) 



For a source region which moves in the time interval 
dt from point A to point B with relativistic velocity v z 
and at an angle ip to the line of sight, the observed differ- 
ence in arrival times for radiation emitted at A and B is 
generally given by dt b s = dt — dt (v z /c) cos ip, thus lead- 
ing to a shortening of the observed time interval. Along 
this line of argument, one may easily derive that the ob- 
server in the model presented here will only perceive a 
strongly shortened period, i.e. the observed period P bs is 
related to the intrinsic period Pk by (cf. also Camenzind 
& Krockenberger 1992, Roland et al. 1994) 



P 



obs 



Pk 



(l + z) / (1-A, cos0(t))di. 



(5) 
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Performing the integration, one immediately arrives at 



Pn 



{l + z)(l-—cosi)P k 



(6) 



From the theoretical point of view, relativistic blazar jets 
are thought to be oriented at a small viewing angle. Cur- 
rent emission models favour an inclination angle i ~ 1 /jb 
(Spada 1999; cf. also Chiaberge et al. 2000) with typical 
bulk Lorentz factors in the range 10 — 15 (e.g. Mannheim 
et al. 1996; Hillas 1999; Spada et al. 1999). For such values 
and by using an observed period of 23 days and a char- 
acteristic outflow velocity of v z /c ~ (1 — I/75) ' 5 , Eq. (^|) 
results in an intrinsic period of = (6 — 14) yrs. 

Combining Eq. (^) and Eq. (^) we may also derive an 
expression for the the centre-of-mass distance 



R = 



Po 



bs 



jl/(3+a) _ J c 



(7) 



2?r (1 + z) /V(3+a) + 1 sini ' 

Given the observed period and the spectral index, Eq. (0) 
only depends on the inclination angle. Accordingly, for 
an observed period of 23 days, a ratio / = 8 and a TeV 
spectral index of a ~ (1.2—1.7) (cf. Aharonianet al. 1999), 
one gets R ~ (2.0 — 3.5) x 10 16 cm using the inclination 
values above. 

By inserting Eq. ([!]) in Eq. (R) , the appropriate binary 
mass ratio is given by 



M 



P. 



1/3 



obs 



(m + M) 2 / 3 (2 7r[l + z]G) 1 /3 S ini 



l/(3+a) 



fl/(3+a) + I 



cos i) 



2/3 



(8) 



For a secondary mass in the range of (10 6 — 10 s ) Mq the 
required primary masses are calculated in Fig. |l| (see the 
curves K (10), -K"(15)) for two different inclination angles 
and / = 8 yielding primary masses of the order of 10 8 Mq. 

2.2. A gravitational constraint on the binary separation 

Observationally, BL Lacs are in general less luminous ra- 
dio sources, showing a lack of strong optical emission 
lines and little signs of cosmological evolution (cf. Bade et 
al. 1998; Cavaliere & Malquori 1999). Celotti et al. (1998) 
have suggested that BL Lac objects correspond to the fi- 
nal evolutionary stage of sources accreting at low radia- 
tive efficiencies (i.e. a dormant black hole system), which 
seems to be supported by HST observation indicating that 
the less luminous AGN stages occur after the original 
quasar has dimmed (Bahcall et al. 1994). Recently, Vil- 
lata & Raiteri (1999) have argued that BL Lacs represent 
advanced and close BBHS with a decreased mass accre- 
tion rate, the binary separation in the case of Mkn 501 
being of the order of that for which gravitational radi- 
ation becomes dominant. Thus, we might set an upper 
limit on the allowed binary masses in Mkn 501 by assum- 
ing that the current separation equals the gravitational 



separation d g , i.e. the position where the gas dynamical 
time scale is balanced by the time scale for gravitational 
energy losses (BBR 80). Gas, which may be constantly 
supplied for example by tidal interaction between galax- 
ies (cf. Heidt 1999) and accreted onto the more massive 
black hole, will cause the binary separation to shrink on 
a time scale i gas ~ M (1M yr -1 / M), with M the accre- 
tion rate (BBR 80). For a simple estimate let us assume 
that during the optical bright QSO phase mass accretion 
occurs at about the Eddington limit. The phase of nuclear 
activity seems to be rather short with a typical duration 
of a few times 10 7 yrs (Haehnelt et al. 1998; Richstone et 
al. 1998). In particular, for a duration of nuclear activity 
of the order of the salpeter lifetime t s — eo~Tc/^TrGm p = 
4.5eo.i x 10 7 yrs and for a primary black hole mass of 
~ 10 8 Mq, gas infall rates of ~ 2MQyr _1 are required 
to sustain the Eddington luminosity, using a canonical 
10% efficiency. On the other hand, assuming a circular 
orbit, the time scale for gravitational radiation is given by 
Tgrav = 6.3 x 10 4 df 6 / (M 8 m 8 [m 8 + M 8 ]) yrs, where the 
distance and the masses are expressed in units of 10 16 cm 
and 10 s Mq, respectively. Thus, by equating the gas dy- 
namical time scale i gas with r grav , the separation at which 
gravitational radiation becomes dominant may be written 
as 



d gi6 = 6.3 M 8 1/2 m\ /4 (m 8 + Ms) 1 / 4 



f M yr- 1 
V M 



1/4 



(9) 



From Eq. ([)]) and Eq. ([!]) we immediately get the relation 



(m + M) 1 / 1 ' 



1.29 x 10 



22 



M 



1/4 



M yr-i 



P 



2/3 



obs 



(2tt[1 + z]) 2 / 3 (1-^ cos?) 2 / 3 



(10) 



This mass dependence is illustrated in Fig. [I] (curves G). 
The respective upper limit is given by the point of inter- 
section with the relevant curve K . For example, applying 
a = 1.2 and using i — 1/10, we have a maximum sec- 
ondary mass m ~ 4 x 10 6 Mq and a corresponding pri- 
mary mass of M ~ 10 8 M (cf. also Table ||). The masses 
shown in Fig. [l] are in a reasonable range for ellipticals. 
Masses of the order of one million solar masses for the 
companion black hole appear to be in agreement with the 
concept that the galaxy swallowed in the merger process 
was a minor spiral galaxy. On the other hand, the host 
galaxy of Mkn 501 seems to belong to these classes of el- 
lipticals which have black holes in the centers of at least 
a few hundred million solar masses. Therefore, the binary 
scenario for Mkn 501 seems not unlikely. 

3. Discussion 

In this paper we have suggested that the periodicity in 
the flaring state observed in Mkn 501 might be caused 
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0.5 1 1.5 



M | 10 a M o ] 

Fig. 1. Required mass dependence for a BBHS in Mkn 
501. The solid [if (10)] and long dashed [if (15)] curve are 
given by the Doppler condition Eq. (||) for inclination an- 
gles i = I/76 with 7f, = 10 and 15, respectively. The curves 
G(l/i) are given by the condition that the current bi- 
nary separation equals the gravitational distance d g , see 
Eq. ( |Io| ) , thus yielding upper limits on the allowed binary 
masses. A TeV spectral index of 1.2 has been applied for 
the calculation. 

by the orbital motion of the jet in a close BBHS. Apply- 
ing a simple toy-model we have shown that the BBHS 
may have a period of ~ (6 — 14) yrs and a centre-of- 
mass distance of ~ (2.0 — 3.5) x 10 16 cm. If one assumes 
that this separation corresponds to the distance at which 
gravitational radiation becomes important, several upper 
limits for the binary masses may be derived. These mass 
ranges, which are shown in Table [j] using an observed pe- 
riod P bs = 23 days seem to be in line with the expecta- 
tions from merger scenarios and the suggestions made by 
Villata & Raiteri (1999). 

The TeV observations indicate that we may have N < 
6 for the number TV of periodic oscillations (cf. Aharonian 
et al. 1999; Catanese & Weekes 1999; Quinn et al. 1999), 
which results in a required propagation length for the 
emitting component of l z = N PkV z ~ 11 — 26 pc. Thus, 
for the projected length at the position of Mkn 501 one 
finds l p ~ 1.4 - 2.1 mas for i = (1/15) - (1/10) rad. Re- 
markably, the jet of Mkn 501 bends dramatically at about 
3 mas from the core (Marscher 1999). Hence, a change in 
the jet parameters might be the reason for the termination 
of the observed periodicity. 

For the proposed model to be valid, the jet has to 
be perfectly collimated with an intrinsic opening angle 
of less then arctan(d/Z z ) ~ 0.05°. Such values are in- 
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i = 1/76 


1/10 


1/15 


m [10 8 M Q ] 


0.040 (0.061) 


0.33 (0.42) 


M [10 s M ] 


0.96 (0.75) 


1.03 (0.91) 


d [10 lb cm] 


2.40 (2.24) 


4.57 (4.54) 


Pk [yrs] 


6.10 


13.7 


Tgrav [10 ' yrs] 


5.46 (4.25) 


5.86 (5.18) 


P P [10 4 yrs] 


0.52 (0.39) 


3.03 (2.93) 



Table 1. Maximum binary masses, separation d, intrinsic 
orbital period Pk, gravitational lifetime r grav and preces- 
sional period P p for inclination angles i, accretion rate 
M = 2Mq/yt and spectral index a = 1.2 (1.7). 



deed expected in scenarios for the formation and col- 
limation of magnetized BL Lac jets (cf. Camenzind & 
Krockenberger 1992; Appl & Camenzind 1993; Schramm 
et al. 1993). At first sight, such a cylindrical jet structure 
seems to be at least ~ 20 times more collimated than the 
jet seen on VLBA maps (cf. Marscher 1999). However, 
there is evidence for an at least two-component jet struc- 
ture in Mkn 501 suggesting an inner spine with transverse 
magnetic field and an envelope with longitudinal magnetic 
field (Aaron 1999; Marscher 1999), the polarization prop- 
erties of the inner spine strongly supporting shocked-jet 
models (cf. Attridge et al. 1999). In fact, our model re- 
quires that the high energy emission originates in a chan- 
nel along the jet axis as in two-fluid models (e.g. Sol et 
al. 1989, Roland et al. 1994), the inner emission probably 
being self-absorbed on the VLBA scale. Recent observa- 
tions of radio jets indeed indicate a confinement of the 
higher energy emission to a well-defined channel within a 
much more extended radio emission (Bahcall et al. 1995; 
Perlman et al. 1999; Swain et al. 1999). The unification 
of BL Lacs and FR I objects may add another piece of 
evidence to such a jet configuration: in order to account 
for the observed spectral properties an at least two-fold 
jet velocity structure seems to be required in which a fast 
spine is surrounded by a slow (but still relativistic) layer 
(Chiaberge et al. 2000). Support for such a possibility is 
positively provided by numerical simulations (cf. Aloy et 
al. 2000; Frank et al. 2000). 

For a TeV flux ratio between minimum and maximum 
of ~ 8, the corresponding shift in the break frequency 
would be given by a factor of ~ (1.5 — 1.6) while the X-ray 
flux ratio becomes ~ (5 — 7) applying an hard X-ray spec- 
tral index a ~ 0.6 — 0.9 (cf. Lamer & Wagner 1998; Pian 
et al. 1998). Such values seem to be consistent with Bep- 
poSAX observations (Pian et al. 1998) and may also be 
recovered, using a broken power law fit, in RXTE observa- 
tions of Mkn 501 (cf. Krawczynski et al. 2000, their Figs.l 
and 2a). Gamma-ray observations carried out by the CAT 
Telescope also reveal a shift of the maximum peak en- 
ergy apparently in accordance with the expectation above 
(Djannati-Atai et al. 1999). Small changes in the maxi- 
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mum electron Lorentz factor or the magnetic field along 
the trajectory of the emission region may further add to 
flux variations. If there is indeed an additional flux contri- 
bution, e.g. low energy emission from the layer, a station- 
ary component comparable to the observed infared-optical 
flux (e.g. Pian et al. 1998; cf. also Kataoka et al. 1999) 
or an additional component responsible for the soft X-ray 
emission (e.g. Lamer & Wagner 1998; Wagner et al. 1999), 
the amplitude of the Doppler modulation may decrease to 
lower frequencies. 

In the simple model presented above, we have not yet 
considered the influence of jet precession due to gravito- 
magnetic and geodetic origin with a period P p ~ 580 x 
d% 2 (Mg+mg) 1 / 2 TOg" 1 M 8 - 1 (l+3M/4m)- 1 yrs (cf. Thorne 
et al. 1986). Since this driving period is much larger than 
the orbital period (cf. Table [j]), a precessional modula- 
tion should be negligible during a few revolutions. Inter- 
estingly, a precessional period of ~ 10 4 yrs agrees with 
the driving frequency found by Conway & Wrobel (1995) 
in order to explain the misalignement of the radio jet in 
Mkn 501 on parsec and kiloparsec scale (see also Villata 
& Raiteri 1999). If the binary hypothesis is correct, the 
observable period should remain similar during different 
outburst phases unless there is a change in the general 
jet properties. For example, an increase in the observed 
period should then be accompanied by a decrease in the 
bulk Lorentz factor or, on larger time scales, by an in- 
crease of the inclination angle due to the jet precession 
(cf. also Eq. g). 
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